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Abstract 
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retardations in the "near zone", where ^ < 1 (rf = oscillation amplitude of the source, 
r = distance from the source). We show that they do not apply in the region ^ ~ 1, and 
prove that gravitational forces are retarded in the immediate vicinity of the source. An 
experiment to measure this retardation is proposed, which may provide the first direct 
experimental observation of propagating gravitational fields. 
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I INTRODUCTION AND SURVEY OF RESULTS 



Between the 1940s and the early 1960s, several researchers proposed the possibility of measuring 
the retardation of gravitational interaction using laboratory experiments [1, 2, 3, 4]. Because of 
the poor technology at that time and contradictory theoretical predictions about the retarda- 
tion, interest in the subject declined. No experiments were performed during this period, and 
after 1962, the topic was never discussed again in the literature. At the 1996 Virgo conference, 
a new experimental set-up was described for the purpose of testing the retardation of gravi- 
tational interaction in the immediate vicinity of a vibrating mass distribution [5]. Discussions 
following the talk indicated that there is still ambiguity about the existence of such a retar- 
dation. In this paper, we will discuss the common arguments which lead to the conclusion 
that gravitational forces generated by oscillating mass distributions are instantaneous in the 
near-field. All of these arguments have the weakness that they are not valid at extremely small 
distances from the oscillating source. We will analyze the gravitational forces in the imme- 
diate vicinity of a specific oscillating mass distribution, for which the dominant part of the 
corresponding linearized Einstein equations can be solved. It shows that there is a possibly 
detectable retardation of the gravitational forces in the immediate vicinity of the source, which 
converges to zero with increasing distance. The dominant length scale for the decay of the re- 
tardation is given by the oscillation amplitude of the mass distribution and can be very small. 
Experimental verification of this local retardation effect would provide the first direct observa- 
tion of propagating gravitational fields. We propose a laboratory experiment which requires a 
highly accurate phase measurement to test our result. Given the state of modern technology, 
this measurement may now be possible. 

II DEFINITION, DISCUSSION AND ANALYSIS OF THE 
MATHEMATICAL MODEL 

We consider an idealized experiment in a laboratory environment, consisting of a mass-balanced, 
mechanical two-body oscillator, which generates disturbances in the space-time metric. The 
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oscillator consists of a point mass to which is connected by way of a spring with negligible 
mass to a second point mass M ^> to. The resulting mass-spring system oscillates with eigen- 
frequency uj about its center of gravity, and does not interact dynamically with additional, 
external masses. The distance between the masses to and M is assumed to be much larger 
than the oscillation amplitude d of the smaller mass to. The oscillation amplitude D = j^d of 
M, which is calculated with classical mechanics, is much smaller than d and leads to negligible 
time-dependent gravitational effects in the immediate vicinity of to, as will be shown in section 
II. 1. We will mainly consider the time-dependent gravitational forces generated by to, and will 
thus refer to to as the "source", or the "source mass". Another mechanical oscillator, tuned 
to resonance to detect these gravitational forces, is placed on the oscillator axis at a distance 
of order 0(d) from to. The purpose of this paper is to show that the gravitational forces gen- 
erated in this extreme near-field laboratory experiment are retarded, which may be verifiable 
with modern technology. Due to many objections against experiments of this kind in the past, 
we will first discuss the common counterarguments. 

The occurrence of a measurable retardation in the proposed system may seem problematic 
for several reasons. One may suspect that the theorem about the non-existence of gravitational 
dipole radiation is violated by the mass-balanced, mechanical two-body oscillator described 
above. It states that isolated mechanical systems cannot generate any dipolar gravitational 
radiation, due to momentum and angular momentum conservation within the system [6] . How- 
ever, this argument does not contradict the existence of retarded gravitational forces which are 
strongly localized around the source. We will refer to the immediate vicinity of to as a region 
around to with a diameter of order 0(d). In section II. 1, we will prove that the time- dependent 
gravitational forces due to the larger mass M can be neglected in the immediate vicinity of the 
source mass to, although the static field due to M may dominate over the static field gener- 
ated by to. The mass to, which is by itself not an isolated system, appears locally to be an 
isolated oscillator. Because the momentum of to is by itself not conserved, to is allowed to 
radiate gravitationally in its immediate vicinity without contradicting the overall conservation 
of momentum. At larger distances from the oscillator, the time-dependent gravitational forces 
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due to M and m have equal strength, and we will show in section II. 2.1 that their radiative 
contributions cancel. Hence, there is no contradiction to overall momentum conservation. 

Near- field set-ups similar to the proposed one are usually analyzed with the "near zone" limit 
used in radiation theory, typically defined as a distance range comparable to the wavelength of 
the radiated signals, which is usually very large. At these large distances from the source, the 
vibrating mass distribution can be very accurately modelled as a series of multipole coefficients, 
which neglect all retardation effects within the source. But in the immediate vicinity of the 
source, i.e. at distances comparable to its vibration amplitude, it is not possible to make 
this simplification. Instead, one must take local retardation effects into account, even if they 
disappear at large distances. 

Another common objection against the existence of observable retardations in the proposed 
set-up arises from an electrodynamic analogy. Retarded electromagnetic potentials may result 
in almost instantaneous electromagnetic fields, due to certain cancellation effects [7] . This phe- 
nomenon occurs in the "near zone" introduced above, and is consistent with our results. We 
will show in section II. 2. 2 that the retardation of the gravitational forces decreases continuously 
to zero with increasing distance from the source, and that in the immediate vicinity of m, the 
retardation has possibly measurable values. 

We will first show that our proposed system can be analyzed with the linearized Einstein 
equations. Let c ps 3- 10 8 [^] denote the vacuum speed of light. We choose a local Minkowskian 
coordinate frame {ct,x,y, z}, in which m oscillates about the origin with amplitude d, angular 
frequency u>, and parallel to the z axis 



z, denotes the unit vector in the z direction, and z(t) is the location of m at time t. We 
assume mass m ~ OfdCT 1 ^]), oscillation amplitude d ~ O(10 _2 [m]), and angular frequency 
id ps O(10 3 [^]) as a set of reasonable parameters. 

To keep track of order- of- magnitude estimates, we introduce dimensionless variables. We 
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define the dimensionless four-vector (r, £) := h (ct,x) and use the oscillation amplitude d to 
represent the length scale. The energy momentum tensor T^ v is related to the dimensionless 
T^ v by =: r ^§-T M1 ,, and the Einstein tensor which has the dimensions [m~ 2 ] introduces 
='■ TpG^v. The Einstein equations, written in terms of dimensionless variables, yield 

8nmG 



in which G ~ 6.67 • 10 11 [ m t 1 is the universal gravitational constant. The number 87r j?T 

*-sec z kg] ° era 

10 -25 estimates the strength of gravitational signals due to m. We bound the mass of M by 
< O(10 2 ), such that all considerations about m also apply to M. The RHS of Eq. ( |2|) is 
very small if || T^ v ||<C 10 25 , and the linearized Einstein equations can be used in this case. 
Consequently, one expands the space-time metric g^ v = r]^ u + around the flat Minkowski 
metric T]^, and introducing 7^ := h^ v — \rj^ v h\, Eq. ( 0) reduces to 

(dl - A f ) 7 ^(r, = - T„(r, , (3) 

by use of the Hilbert gauge condition 7^ = [8]. A|*is the Laplace operator in the dimension- 
less spatial coordinates £. Eq. ( ^) is solved by the superposition of a particular solution with a 
homogenous solution. The full particular solution includes the influences of both M and m. We 
impose the boundary condition 7^ — > in the limit |£| — > 00. Since the full particular solution 
is asymptotically zero, as will be shown in section II. 2, this requires the homogenous solution, 
given by a superposition of plane waves with constant amplitudes, to be zero. Therefore, 7 M „ is 
the particular solution of Eq. ( |3|). Due to linearity, the contributions from M and m to 7^ and 
T^ v are additive and can be discussed separately, thus 7^ = 7^ + 7^ and T M „ = Tffl +Tf£f\ 
We will first analyze the influences due to m, determined by 7^ and . We obtain 



7^(r,0 = - J dr>d 3 e ^_ p 1 5{r' - r + \£ - ?\) 0(r - r>) , (4) 



in which 5 is the Dirac-delta distribution and is the Heaviside function. The perturbation 
parameter of our analysis is the ratio (3 = — ~ O(10 -8 ) between the maximal velocity ud of 
the oscillating mass m and the speed of light c. Our analysis will show that the retardation is 



an effect of the order 0(/3). Therefore, a relative accuracy of the order 0((3 2 ) is sufficient for 
our calculations. Tqo can be decomposed into 

where T g is of order 0(/3*). Note that this decomposition does not correspond to the full 
Taylor expansion with respect to (3. Tq™' ^ is the only contribution of order 0(1 = (3°) to 
T MI/ , and is determined by the dimensionless rest energy density p(t',£') of m, which satisfies 
/ d(,' p(V , £') = 1. In all physical systems, p(t', £') can be assumed to be a sharply-peaked, 
smooth function with sup^, p(t',£') <C O(10 25 ), such that Eq. ( [|) can be used. For math- 
ematical convenience, we replace p(t',£') by a Dirac delta distribution 6^ 3 \ and consequently 
obtain 

T { o mi0) (r',e / )=^ (3) (e'-^in/3r'). (5) 

Note that the space and time dependent terms on the RHS of Eq. (5) cannot be separated, 
as opposed to the expression p(£') expficur'] commonly found in the literature which is suitable 
for multipole decomposition. Tq™' 2 ' ) is determined by the sum of the kinetic energy of m and 
the potential energy of the spring, and is of order 0{j3 2 ). The components , i = 1,2, 3, are 
given by 

Tt\r',i') =Pz t cos/V - zsin/V) , 

and are of order 0{j3). The components i,j = 1, 2, 3, given by 

T^\t',(') = (3 2 Zi % cos 2 f3r' 5 {3 \e - zsin/V) , 

will only contribute to error terms of the order 0(/3 2 ). 

We evaluate 7^ (r, £) at the point ( = on the positive z axis, such that £3 = ^ > 1. 
Inserting Eq. (5) into Eq. (4), we first calculate 7oo 1 '°^(' r > £3^)5 defined by 700 — 7ocT'°' ) +7ocT' 2 ' ) - 
The 5-integration in ^ and £2 can be readily carried out, yielding 
(m,o)/_ t ^ 4mG f 6(63 -sin/V) 



7^' Uj (r, ^) = - y VSJ _ ^ ' 5(r' + $3 - ^ - T ) 0(r - r') 



c 2 d J £3 — sin f3r' 
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Further reduction is possible using the formula 

fdr'5(f(T'))g(r')= £ -i- g{6) , 

ee{<?|/(0)=o} ^ ^ ' 

in which we identify the zeros of the argument in the delta distribution as the roots of 9 in the 
equation 

0(r,e 3 )=sin^(r,6)+r-6 • (7) 

Note that #(t, £3) is a function of the combination r — £3. The derivative of the argument of 
the delta in Eq. (6) yields 1 — (3 cos /38(r, £ 3 ) at 0(r, £ 3 ). We obtain for 7oo 1 "°' ) (' r ? £3^) 

4mG 1 1 



(m,0) / c a\ 
7oo fa 53*0 



c 2 d l-/3co S/ 30fa£ 3 ) £ 3 -sin/30fa£ 3 ) 

In Eq. ( H), /3 is not only a small perturbation parameter, but also a " dimensionless angular 
frequency". In contrast to j3, the arguments [39(t, £ 3 ) of the trigonometric functions need not 
be small for arbitrary r or £ 3 , and cannot be used as perturbation parameters. Since Tq™' 2 '* 
is of order 0([3 2 ), it follows that 7q™' 2 ' ) = 0((3 2 ) is an error term, and that 7o[T' ) fa£^) = 
1oo'°\ T i&) + 0((3 2 ). Hence it is sufficient to analyze 7o™'°' ) to the order of 0(1 = j3°) and 
0(f3), to determine 7qq with a relative error of order 0{(3 2 ). The solutions for 7Q™' ) fa £ 3 i:), 
i — 1, 2, 3, can be derived in a similar manner, yielding 



(m) / ^ ~\ 

7oi fa 6^) 



4mG Zi /?cos/?#(t,£ 3 ) 1 



c 2 d 1-/3 cos /3flfa£ 3 ) £ 3 -sin/?0fa£ 3 



The solutions for 7 t -]: fa, £32) , i, j = 1,2,3, which are given by 



(m) 4mG zt% /3 2 cos 2 /3#(r,£ 3 ) 1 



c 2 d 1-/3 cos /?0fa,£ 3 ) £ 3 -sin/?0fa,£ 3 ) ' 
contribute to error terms of order 0{(3 2 ). The contributions 7«^fa, £3^) of the mass M can 
be obtained by substituting d-*D,r^R, m^M,(3-^ # ~ > fa>£) — > fjfa, £), 

£3 = 2 ~ * % i n 7u™K r > £3^)5 where .D = is the oscillation amplitude of M, and i? is its 
average distance to the point of measurement. This is because both masses oscillate sinusoidally 
with the same frequency u and opposite phases along the z axis. 
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It follows that 9 r 7o7'' (t,£ 3 z) is of order 0([3 2 ) for both M and m. We will show below 
that the only Levi-Civita connection coefficients needed for our analysis are 

r oo( r >^) = -^(2V,oO"»&2) - 

= 1 7oo, M ( r , 6^) - 7o m ,o(t, £ 3 z) + \ Ik At, &£) . (9) 

Thus, Tq = ^7qo^ + 0((3 2 ) for i = 1,2,3. Here, we have used that = 7^ — |^7a and 
^pa = ^V^ u (h P u,a + h au p — h aP;l/ ), which include the influences of both M and m. 

Given the metric g pu = rj^ + h^, we intend to study the motion of a test mass point rh Cm 
in the immediate vicinity of m, that is part of another mechanical two-body oscillator on the 
z axis. The dimensionless space-time coordinates of rh are scaled with d, yielding (£°(t), |*(t)). 
We are only interested in oscillatory solutions with small amplitudes of order d ■ 0(13), 
which are induced by the oscillations of the space-time metric. Therefore, we require that the 
dimensionless velocity components £* , i = 1,2,3, of rh, are of order 0(f3 2 ). The subscript ",o" 
denotes the derivative with respect to r. The proper time of the test mass satisfies £°(t) t 
and £° = 0(1), since the deviations from flat Minkowskian space-time induced by M and m 
are very small. The action which determines the dynamics of the mass rh in our coordinate 
system is given by 

S = rhcd£ dr{-^^ - U(C)} , 

where £/(£*) is the dimensionless potential energy due to the spring. For the moment, we 
assume that the second mass of the pick-up system is "infinitely" larger than rh, and static. 
The space-time trajectory of the test mass rh is determined by the Euler-Lagrange equation 

£,00 + ^pa^fi +J j _ q 

where are linearized Levi-Civita connection coefficients. The dominant terms in the sum 

r^fo£,o arise from r oo£°oCo • A11 terms I^^o which involve €fo, 3 = 1,2,3, are at least 
a factor of order 0(f3 2 ) smaller, and can be neglected. From the equation of motion for £° 
follows that + r[] = with a relative error of order 0((3 2 ). Using the fact that r[j is a 
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total derivative with respect to r, and that || 7^ ||<C 0(/3 2 ), one obtains £° = 1 + 0{f3 2 ) and 
sjg^viflifi = 1 + 0(/9 2 ). Consequently, the motion of m is determined by 

6,00 + ^ + ^7^ = 0, (10) 

with a relative error of the order 0{j3 2 ). The contributions of M and m to the "gravitational 
force" 1 7oo|j are additive, and have the same functional form. This weak-curvature limit is 
usually referred to as the post-Newtonian approximation, or the Newtonian limit. The pick-up 
oscillator is typically driven at resonance, but has a finite Q factor. In all realistic systems, the 
condition that the velocity is of order « 0{(3 2 ) is satisfied, because the gravitational excitation 
of the pick-up system is very weak, and Eq. ( 0) can be used. 

The Levi-Civita connection coefficients Eq. ( |9|) depend on our specific choice of coordinates, 
or gauge. It is not a priori clear that these coefficients are the correct expression for the physical 
gravitational force, since they could be influenced by meaningless coordinate effects. The 
physically measurable force is determined by the change of these coefficients in the absence of the 
source, i.e. for T^ v = 0, calculated in the same gauge, and using the same boundary conditions. 
In the Hilbert gauge, T^ v = results in the vacuum field equations (d 2 — A^j 7 Al „(r, |) = 0. 
The solution is a superposition of plane waves with constant amplitudes, and must be zero due 
to the boundary condition 7 Mi/ — > for |£| — ► 00. Consequently, the corresponding Levi-Civita 
connection coefficients are zero, which shows that the terms Tq in Eq. ( % = 1, 2, 3, correctly 
describe the physical force. 

In the Newtonian limit, times (8) is interpreted as the classical gravitational potential 

, „s _ mG 1 1 . . 

(T '^° = ~~ l-/?cos/?0(T,f 3 ) f3-sin/?0(T,f 3 ) ( j 

generated by m, which acts on the external mass point m by a gravitational force that is 

proportional to the gradient of V, according to Eq. ( |10D . The physical meaning of the solution 

Eq. ( [TT]) becomes more apparent if one expands the first fraction in a geometric series, and 

the second fraction in a Fourier series 

mG 



V{t,rz)= -— \l + (3 cos f39(- t, -) + 0((3 2 )} 
r I da ) 



E (;) 2 1 - [ l + c n ] smn(3 6( C -t } r -) + 

n>2, n odd V / 

E (-) 2 1 - [ l + d n ] cosnpO^t, r -) \ , (12) 

n>l,neven \ r / GS GS J 

in which all coefficients c n and d n are of order 0((f,) 2 ). The leading terms of order 0(1 = j3°) can 
be obtained by setting (3 equal to zero in the first bracket. The time- independent contribution 
^{1 — {^j } _1 / 2 in the Fourier expansion is the static potential of a resting mass distribution 
centered at the origin, and converges to the classical Newtonian potential — for large r. The 
next, first harmonic term has the asymptotic form 2^(1 + 0((^) 2 )) sin/3#(^t, ~). For large r, 
it describes a propagating, sinusoidal oscillation proportional to \, which decays more rapidly 
than the usual radiative - solutions. All following terms are higher harmonics that decay 
even faster with increasing distance. The contributions of 0(f3) exhibit the same asymptotic 
behaviour with respect to -. 

II. 1 GRAVITATIONAL EFFECTS OF THE MASS M 

In this section, we will discuss the gravitational effects of the mass M in the immediate vicinity 
of m. The average distance R between M and the point of measurement, lying in the immediate 
vicinity of m, is much larger than d, and also much larger than r. Due to M > m, the oscillation 
amplitude D = j^d of M, which is calculated from classical mechanics, is much smaller than 
d. Both masses oscillate at frequency uj along the z axis, and the gravitational potential due 
to M at £ = 2% can De obtained by inserting d — >• D, r — > R, m — >• M, [3 — > —^(3, 9 — > j^O, 
(t, £) — > |j(t, £), £ 3 = ^ — >■ ^ in Eq. ( [L2D. The static term is approximately — and 
can dominate over the static field stemming from m. However, one obtains approximately 
^Ji~Ji = !n fitr(ji) 2 f° r the amplitude of the first time-dependent term in Eq. ( [Tjp, which is 
a factor of (-^) 2 smaller than the corresponding term due to m. The amplitude of the first 
harmonic term due to M has the opposite sign as compared to corresponding term due to m, 
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because M and m oscillate with opposite phases. 

In addition, the amplitudes of the higher harmonic terms, approximately given by ^p(^) n = 
r ( ) n (^) n+1 (^)" _1 ' are a f ac t° r °f (■^)™ +1 (^) ra_1 smaller than the corresponding terms due 
to m. One can therefore neglect the time- dependent gravitational influences due to M if is 
small, which is the case in the immediate vicinity of m. At larger distances from m, the number 
approaches 1, and all radiative effects cancel, which is required by momentum conservation. 

II. 2 RETARDATION OF GRAVITATIONAL FORCES 

It will first be shown that the forces derived from the gravitational potential Eq. ( [IT]) are 
not retarded in the region known as the "near zone" in the literature, which is in accordance 
to [7]. Then, it will be proved that these gravitational forces are retarded in the immediate 
vicinity of m. As discussed in the derivation of Eq. ( (TJJ), we place a small test mass fh which 
is part of a resonant mechanical oscillator at position £ = r z on the z axis to detect the 
gravitational forces generated by the source m. From F(t,^z) = —md r V(T,^z) it follows 
that the gravitational force acting on fh is 

F(r, ~J) = - ^ ( ( - R * M - . * r .X + 0(A . (13) 
d d 2 y\l - (3 cos (3e{r, 2 ) ^-sm^(r, 2 ) J J 

F(t, 2%) points along the z axis, whereas the other force components vanish due to symmetry 
reasons. Eq. ( |13|) is accurate to orders 0(1 = /3°) and 0(f3), but not to the order 0(f3 2 ), since 
Tq™' 2 ' 1 has not been considered in 7o™'°' 1 . For our analysis it suffices that the results derived 
from Eq. ( 0) are conclusive to 0(1 = f3°) and 0{(3). To obtain Eq. ( [13]), we have used the 
fact that the derivative of 9(t, ^) with respect to r is given by 

r 1 1 

d r 9(r,-) = -- - _ - - ^- g (14) 

which follows from 

r v 1 

dr 0(t, -) = d r sin /3 6(t, -) - - 

a d d 

= (3 cospe(r^) drOirJ-) - - . (15) 
a d d 
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Note that the leading term in F(r, ~) is proportional to V(t,^z) 2 . The expression for the 
gravitational force generated by M is similar to Eq. ( |l^). It has a static component which is 
much stronger than the static term due to m. However, the time- dependent gravitational forces 
generated by M are negligible in the immediate vicinity of m, because - ~ 1 and therefore 
^ < 1. At large distances from m, M and m produce time- dependent gravitational forces of 
the same strength, since - — > and 4j — > 1. 

° ' r R 

II. 2.1 GRAVITATIONAL FORCES IN THE "NEAR ZONE" 

The spherical region {r = \x\] (3 <C - <C 1} is commonly defined as the "near zone" in the 
literature. For this region, the approximation @{t,^) ~ T — \ can be used, which is familiar 
from radiation theory. In addition, we have used the identity = 1 — y + 0(\y\ 2 ), with \y\ < 1, 
to simplify the fractions in V(t, '^z), and to reduce the potential to 

V(t,-J) = (l+/3cos/3(r \ {-+(-) sin Hr -) } + O{0 1 ) . (.1.0) 



<i d I gT 

The trigonometric functions can be approximated by their Taylor expansions to first order in 
the variable ^ C 1. Keeping terms only to the orders 0{(3) and O(^), one obtains 

F(tJ-z) = -H^l {l + 2/3 cos M |l + 2-sin/?r-2/? cos/?rl + 0((3 2 ) 
a r l [ r J 

= — ll + 2- sm /3t + 4/3 - cos f3r em fir > + 0(/r) (17) 

f I 7" 7" I 

for the gravitational force acting on the test mass m. The term of order 0(/3f) <S is a 

negligible contribution to the second harmonics. Hence the contributions of order O(0) cancel. 
The time-dependent term of order 0(1 = /3°) in the force acts instantaneously, thus there is no 
measurable retardation in the "near zone" where - < 1. The electromagnetic version of this 
result can be found in [7]. Historically, this result has convinced the physics community that 
there exists no region at all where any retardation can be observed. 

The additional constraint ^ -C 1 allows the influence of the second mass M to be neglected. 
As r increases, this condition can not be satisfied, and Jj — > 1. The complete gravitational 
potential in the "near zone" limit, which includes the contributions of both m and M is given 
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by 

j, i r ^ mG L Mr] mGd f . Q fr\ 2 . \ 

which follows from the discussion in section II. 1. The time-dependent part of this solution 
vanishes in the limit — > 1, and there remains a static solution which is dominated by the 
larger mass M. Note that ( m+ ^ G [ s the rotation symmetric limit solution for r — > oo, which 
corresponds to the asymptotically vanishing particular solution of Eq. ( |]). We conclude that 
our result is consistent with [7] , and that it is in accordance with the conservation of momentum. 

II. 2. 2 GRAVITATIONAL FORCES IN THE IMMEDIATE VICINITY OF THE 
SOURCE 

The immediate vicinity of the source is characterized by - ~ 1 , in contrast to - <C 1 which is a 
necessary condition for radiation theory. Consequently, radiation theory breaks down. At this 
distance range, | < 1 holds, and time-dependent influences of M are negligible, as shown in 
section II. 1. As opposed to the situation in the "near zone", it is not possible to approximate 
functions of - by the leading terms in their Taylor series. For this reason, we will have to 
explicitly calculate the dominant Fourier components of the lowest harmonics in Eq. ( |I3|). To 
order 0({3), we find 

The condition f3 - which has defined the "near zone" is also valid for the immediate neigh- 
borhood of the oscillating mass point. Therefore, it is still possible to expand the trigonometric 
functions with respect to (3 5 ~ 0(@). From the defining equation Eq. ( ^) for 9(r, 5), it can 
be seen that (3 9(t, ^) is a function of the variable (3 5. From Eq. ( |14|), one finds 

P0{t, T -) = (3 9(t,0) - P T - + 0(/3 2 ) (19) 

for the leading terms of the Taylor series in (3 ^. We expand the trigonometric functions in Eq. 
( 0) with respect to (3 9(t, ^) « (3 9(t, 0) — (3 5, and perform the Taylor series to first order 
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in (3 in the second bracket. Consequently, the RHS of Eq. ( [L8|) yields 



;i + 2/3cos/30(r, 0)) 



1 ^ (3 cos /30(r,O) 

v l- Jsin^(r,0) " (l_^ si n^(r,0) N2 

The expansion of this product to order 0(/3) gives 

mmG I 1 J/3 sin2/30(r,O) 



+ 0(0 



2 i 



+ 0((3 2 



(l- Jsin^(r,0)) (l-Jsin^(r,0); 
Calculation of the Fourier components of the first harmonics to order O(0) results in 

777777 (J I [ /7 /"/ ] 

-P r h + 2 - sin /Mr, 0) - 2/3$(-) cos/30(r, 0) + /iA + 0(/3 2 ) , (20) 

r ^ _ (<l) 2 V I r r J 

where the abbreviation /i./i. stands for "higher harmonics". The function $ (x) is defined by 

*(*) = ! ((l-xf - 1 + (21) 

for x = J G [0, 1]. $ (x) is a monotonically increasing function on the unit interval with values 
$(0) = and $(1) = 1 at the boundaries of its domain, see Fig. 1. We have only calculated 
the first harmonic terms of maximal order 0(/3) in Eq. ( |2"0| ), because we assume that the 
measurement in our proposed experiment is only sensitive to the first harmonics. Using the 
relationship 

P0{t, = P9(t, 0) - + 0(/3 2 ) , (22) 

which is a generalization of Eq. ( |19D, we finally obtain 

F(t,-z) = -r\l + 2 -sin/30 r, $ - - + /U. + 0(/3 2 ) . (23) 

(!"(?)) ' 
This result shows that the first harmonic term in the gravitational force of order 0(0) can be 

absorbed into the contribution of order 0(1 = (3°). Moreover, it is obvious from Eq. ( |7D that 

9(r, is a function of the combination r — $(-)^, which shows that the gravitational force 

is retarded. For two points 7*1 £ and r 2 z on the oscillator axis, the relative phase difference A0 

of the gravitational force is given by 

r\ a r 2 a 

= - $ - n - $ - r 2 , 24 

c r 1 r 2 
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and is of order 0(/3). At very small distances r from the source, where - — > 1 and $(-) — > 1, 
the value of the phase shift approximates the result which would occur for a wave that travels at 
speed of light. At larger distances in the "near zone", we have - — > and S®( ) ~ > 0- Hence 
there is no measurable phase shift, and the gravitational forces seem to act instantaneously, as 
was shown above. In the intermediate region where - ~ 1, the phase shift can be calculated 
from Eq. ( ^4|). Note that in the limit - — > 0, ( |23[) converges to ( |I7D , as it should be. The 
absolute phase shift 0(r) := ^$(^)r converges rapidly to zero for r — > oo, see Fig. 2. 

We conclude that there is a possibly measurable retardation of the gravitational forces in the 
immediate vicinity of the oscillating source in our proposed system. A discussion of extended 
gravitating sources with the linearized Einstein equations would be a straightforward general- 
ization of the presented analysis. Due to the superposition principle, expressions similar to Eq. 
( |23|) would then have to be integrated over the spatial spread of the oscillating continuum. 
However, we have assumed cylindrical symmetry in our mathematical model, and thus could 
find a solution for the dominating gravitational force term. For systems without cylindrical 
symmetry, the analysis will be much more difficult. From the results in this section, we will 
now discuss how they could be experimentally verified. 

Ill EXPERIMENTAL POSSIBILITIES 

Eq. ( |23|) suggests that it may be possible to experimentally verify the retardation of gravita- 
tional forces in the immediate vicinity of an oscillating mass distribution. This would require 
both the amplitude and phase detection of gravitationally transmitted vibrations. With the 
previous mathematical model in mind, we propose the following experiment: 

A mass m within a mass-balanced apparatus is oscillated with a very stable phase and 
frequency to, as described in the mathematical model. If a mechanical resonator with eigenfre- 
quency uj, mass m, and quality factor Q is placed at position rz, it will be driven into resonance 
by the first harmonic term in the gravitational force. Given suitably chosen parameters, the 
amplitude of the sinusoidally vibrating resonator, typically of the order O(10 _9 [m]), can be 
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detected with modern technology [5]. 

We propose to measure the retardation of the gravitational force in the following way. Place 
the mechanical resonator at the location T\Z on the oscillator axis, such that — ~ 1. From 
the observed gravitationally induced sinusoidal vibration, the phase difference A<pi between 
the phase of the resonator and the phase of the vibrating source mass can be measured. If 
the mechanical resonator is placed at another location r 2 z on the oscillator axis, a different 
result A02 will be observed. By subtracting A0i from A02, all internal influences from the 
experimental measuring instruments can be eliminated, resulting in the relative phase shift 
A0i — A02- The value of c in Eq. ( [24]) can then be calculated from the measured relative 
phase shift A0 2 — A0 1; which would be the speed of light according to General Relativity. 
The expected values for the relative phase shift are shown in Fig. 3. Since linearized General 
Relativity is formally a special relativistic field theory, measurement of c with the proposed 
system would provide a new test of its Lorentz covariance. Moreover, our proposed laboratory 
experiment would provide the first direct observation of propagating gravitational fields. In 
addition, the scalar component of order 0(1 = j3°) in Eq. ( |3|) also determines the Coulomb 
potential of an oscillating point charge after a substitution of parameters. Therefore, our 
solution Eq. ( ^) also applies to the retarded Coulomb force generated by an oscillating point 
charge. Measurement of the retardation in the immediate vicinity of the charge would thus be 
another interesting experimental possibility. 

Since the 1960s, several researchers have detected the amplitudes of gravitationally induced 
vibrations in laboratory experiments [5, 9, 10, 11]. It may be possible that these experiments can 
be analyzed in a similar way as the model in this paper, which would provide information about 
the magnitude of the retardation, if present. The detection of the retardation of gravitational 
interaction in laboratory experiments has not been discussed much recently because the usual 
"near zone" calculations predict instantaneous interactions, as was shown in section II. However, 
this result does not apply in the immediate vicinity of the oscillating source, as was also shown 
in section II. Note that our experiment is not related to the gravitational waves of helicity 
two, which are far more difficult to detect, since they are of order 0{(3 2 ). Unfortunately, 
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several technical problems must be solved before measurements to test these retardations can be 
realized. In 1958, Q. Kerns demonstrated an analog electronic circuit capable of measuring 10 -6 
degree phase shifts [4] . We have made several modifications to Kerns' original design, and we are 
currently able to measure 10~ 7 degree relative phase shifts over a 30 second measurement period 
[13]. Currently, the technology is thermally limited, but thermal control of the system and 
application of noise reduction techniques may soon enable us tho increase the phase sensitivity 
to the 1CT 9 degree phase accuracy required, to measure the retardation to 1% [5]. It is not 
yet known if the phase stability of the mechanical systems used for the experiments can be 
maintained to this accuracy over the measurement time. 

IV CONCLUSIONS 

Using the linearized Einstein equations, we have proved that the gravitational force generated by 
a specific oscillating mass distribution is retarded in its immediate vicinity. We have proposed 
an experiment to measure this retardation which may enable the first direct experimental obser- 
vation of propagating gravitational fields. The amplitudes of gravitationally induced vibrations 
in systems similar to ours have been successfully detected since the 1960s. The detection of the 
retardation requires a very sensitive phase measurement, which may now be possible given the 
state of modern technology. 
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V Figure captions 

Fig. 1: The function <&(x) versus x, where x — f. 

Fig. 2: The absolute phase shift 0(r) for r between d and l[m], at fixed parameter values 
d = 0.01 [m], oo = 251. 3[r ad /sec], and c = 3 • 10 8 [m/sec]. 

Fig. 3: The change in phase shift Acf) for r 2 between r\ and 0.1[m], and for fixed parame- 
ter values T\ = 0.03[m], d = 0.01[m], uo = 251. 3[r ad /sec], c = 3 • 10 8 [m/sec]. Note that the 
value r 1 = 0.03 [m] corresponds to the smallest possible separation between the centers of the 
source mass and the detector mass in our experiment, due to the finite size of the physical 
systems being used. 
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